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Extensive electronic structure calculations are performed to obtain the stable geometries of metals 
like Al, Ga and In on the Si(OOl) surface at 0.5 ML and 1 ML coverages. Our results coupled with 
previous theoretical findings explain the recent experimental data in a comprehensive fashion. At 
low coverages, as shown by previous works, 'Para' dimers give the lowest energy structure. With 
increasing coverage beyond 0.5 ML, 'Ortho' dimers become part of low energy configurations leading 
toward a 'Para' to 'Ortho' transition at 1 ML coverage. For In mixed staggered dimers ('Ortho' 
and 'Para') give the lowest energy configuration. For Ga, mixed dimers are non-staggered, while 
for Al 'Para' to 'Ortho' transition of dimers is complete. Thus at intermediate coverages between 
0.5 and 1 ML, the 'Ortho' and 'Para' dimers may coexist on the surface. Consequently, this may 
be an explanation of the fact that the experimental observations can be successfully interpreted 
using either orientation. A supported zigzag structure at 0.5 ML, which resembles (CH)x, does not 
undergo a dimerization transition, and hence stays semi-metallic. Also, unlike (CH)x the soliton 
formation is ruled out for this structure. 

PACS numbers: 73.63.-b, 73.90.+f, 68.90.+g 



I. INTRODUCTION 

The study of metals on semiconductors dates back to 
the nineteenth century and has seen a vigorous recent re- 
vival due to tremendous interest in Nanotechnology. Our 
ability to manipulate atoms, placing them at will on dif- 
ferent surface sites to create exotic artificial structures, 
has led to further investigations^ of electronic and trans- 
port properties of free and supported nanowires. The 
placement of trivalent atoms (Al, Ga and In) on Si(OOl) 
can lead to the formation of low-dimensional structures, 
exhibiting significant new electronic properties. 

One can easily compute the electronic properties of 
free standing nanowires. For this, one must first com- 
pute the total energy and determine the possible stable 
structures. Such calculations have indeed been carried 
out for nanowire32iiiiSiSiL& consisting of a wide variety 
of atoms, e.g., K, Al, Cu, Ni, Au and Si. A general find- 
ing is that a zigzag structure in the form of an equilateral 
triangle is the most stable2*iL&. This can be understood 
as arising primarily due to the maximization of coordina- 
tion number for each atom in a quasi ID structure. An- 
other structure which also has a local minimum, but not 
terribly stable, is a wide angle isosceles triangle which 
somehow is reminiscent of the bulk environment. For 
example. Si which is a four fold coordinated in the bulk 
(tetrahedral angle 109°) shows^ a local minimum at an 
angle of ~ 117°. In general, free standing nanowires tend 
to be metallic (have bands crossing the Fermi level) but 
these nanowires in practice are to be supported. Silicon 
is the most widely used substrate for practical applica- 
tions and the low index surfaces, Si(OOl) is the surface of 
choice. With the downward spiral toward nano devices, 
it is desirable to study the metallic properties at the low- 
est possible coverages. It is in this context that the study 
of metals like Al, Ga and In at submonolayer coverages 



on Si(OOl) take on the added importance. The interac- 
tion of metal nanowires with substrate can significantly 
alter the electronic properties. 

Low-energy-electron-diffraction (LEED), Auger elec- 
tron spectroscopy (AES) and Scanning tunneling mi- 
croscopy (STM) studies have provided much informa- 
tions about the interaction of Al, Ga and In overlayers on 

the Si(OOl) surface at different coverages and at different 
temperaturea9iiodiii2,i3,i4,i5,i6,i7,i8,i9,20^ 

Ide et al.^ observed 2x2, 2x3, 4x5, and 1 x 7 structures 
depending on the coverages less than 1 monolayer (ML) 
of Al on Si(OOl) and the substrate temperature in their 
LEED and AES experiments. However, the 2x2 and 2x3 
phases were observed around 0.5 and 1/3 ML coverages 
respectively at low temperatures. 

Sakamoto and KawanamiiSi performed Reflection high 
energy electron diffraction experiments and established 
the existence of phases with various symmetries (2x3, 
2 x5, 2x2 and 8x1) for Ga coverages less than 1 ML 
and at temperatures ^ 350° C. Later, Bourguignon et 
aL'^'^ examined the evolution of IML of Ga on Si(OOl) 
with LEED and observed all the above mentioned phases 
and in addition they also observed 1x2 phase at 1 ML 
coverage. Nogami et a L^^d^ did STM studies and ob- 
served 2x2 at and below 0.5 ML coverage of Al, Ga and 
In on Si(OOl). Knall et al>i^ performed experiment for 
In adsorption on Si(OOl) using RHED, LEED, AES and 
STM and observed transition from 2x2 to 2x1 struc- 
ture of In as coverage increased. All the experimental 

resulta'^d°i"d^d3,14,15,16,17,18,19,20 pojj^^ ti^g f^^^ ^^^^ 

at low temperatures Al, Ga and In form 2x2 structure 
on the Si(OOl) surface at 0.5 ML coverage. 

At higher metal coverage, Nagomi et ali^ observed 
3D island formation of Al on Si(OOl) around 1 ML cover- 
age while ordered structures were observed for Ga and 
In on Si(OOl) at 1 ML coverage without any island 
formationiiii^. In particular, a 1x2 structure was ob- 
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served for Ga, and 2x2, 2x1 structures were observed 
for In on Si(OOl). Since Ga does not form islands at 
IML coverage and it has been argued that Al forming is- 
lands at 1 ML coverage is unusual, we present a complete 
theoretical study for Al/Si(001) at 1 ML coverage. 

Note that several experimental results were 
explainec&iiiiSiiiiSiiSiSi in terms of metal dimers 
oriented parallel to the Si dimer rows i.e. metal dimers 
were assumed to be directed perpendicular to the Si 
dimers in the reconstructed surface. This orientation 
of metal dimers is named as 'Ortho' dimers. Electronic 
structural calculations were done by BatraSi for Al 
on Si(OOl) at 0.5 ML coverage and he explained the 
experimentally observed patterns in terms of these 
'Ortho' dimers. Later, Northrup et al.^^ did theoretical 
calculations on the same system and showed that the 
orientation of metal dimers parallel to the Si dimers 
was most favorable. These dimers orientated parallel to 
the Si dimers are called as 'Para' dimers. A series of 
experimentsiiiSiiS performed recently for the structural 
configuration of metal dimers on the Si(OOl) surface 
indeed confirm the findings of Northrup et ali^S at 0.5 
ML coverage. To the best of our knowledge, there is 
hardly any detailed theoretical study at higher coverages 
(1 ML) for these systems. 

In this paper, we put significant effort to perform zero 
temperature electronic structural calculations for both 
0.5 ML and 1 ML coverages of Al, Ga and In on the 
Si(OOl) surface and compare our results with the exper- 
imental and theoretical findings where available. We at- 
tempt to put the work in perspective and fill the holes 
by presenting new results. In particular, we find that 
our calculations at 0.5 ML give results in accord with 
Northrup's predictions and also with the recent exper- 
imental results obtained by Sakama et aliLiSiiS. More 
significantly, we obtain an orientational transition ('Para' 
to 'Ortho') for Al dimers as we go from 0.5 ML to 1 ML 
coverage. Our calculations also reveal that the recon- 
struction of Si(OOl) is completely lifted at 1 ML coverage 
of Al. We find that In and Ga interacts differently with 
the surface; for these metals the reconstruction is not 
lifted at 1 ML coverage. One finds the staggered mixed 
('Ortho' and 'Para') dimers to be the lowest energy struc- 
ture for In while the mixed dimers are non-staggered for 
Ga. The important conclusion is that the 'Ortho' dimers, 
which are energetically unfavorable at 0.5 ML coverage, 
become viable at higher coverage facilitating 'Para' to 
'Ortho' transition as a function of coverage. Also, a pre- 
viously reported semi-metallic phase of supported zigzag 
Al nanowire is found to be stable against a dimerizing 
Peierls transition unlike (GH)x. The soliton22r formation 
in the supported nanowire is argued to be improbable. 

The paper is organized as follows. The calculational 
parameters are given in sec. II and the results and dis- 
cussions are presented in sec. Ill followed by a summary 
of our findings in sec. IV. 



II. METHOD 

First principle total energy calculations were carried 
out within density functional theory at zero temper- 
ature using the VASP code^^. The wave functions 
are expressed by plane waves with the cutoff energy 
fc-l-Gp < 250 eV. The Brillouin Zone (BZ) integrations 
are performed by using the Monkhorst-Pack scheme with 
4x4x1 fc-point meshes for 2x2 primitive cells. Ions are 
represented by ultra-soft Vanderbilt type pseudopoten- 
tials and results for fully relaxed atomic structures are 
obtained using the generalized gradient approximation 
(GGA). The preconditioned conjugate gradient method 
is used for the wave function optimization and the con- 
jugate gradient method for ionic relaxation. 

The Si(OOl) surface is represented by a repeated slab 
geometry. Each slab contains five Si atomic planes with 
Hydrogen atoms passivating the Si atoms at the bottom 
of the slab. Consecutive slabs are separated by a vacuum 
space of 9 A. The Si atoms on the top four layer of the 
slab are allowed to relax while those in the bottom layer 
of the slab and the passivating Hydrogen atoms are kept 
fixed to simulate the bulk like termination. The conver- 
gence with respect to the energy cutoff and the number 
of k points for similar structures has been established24. 

III. RESULTS AND DISCUSSIONS 

We have performed total energy calculations using 
the above parameters for both the ideal and the recon- 
structed Si (001) surfaces. For the ideal case, the Si 
atoms on the surface are arranged in a square pattern of 
side 3.84 A. For the reconstructed surface the Si atoms 
rearrange in such a way that the top layer Si atoms form 
dimer rows along the y [110] direction but the atoms 
move along the x [110] direction to form dimers. Here 
we consider 2x2 supercell for our calculations. At the 
0.5 ML coverage of Al, Ga and In we start only with the 
reconstructed surface and at 1 ML coverage we start with 
both the ideal and reconstructed surfaces. 



A. 0.5 ML of coverage 

Here we consider various structural arrangements of 
metal atoms on the Si(OOl) surface. The 2x2 primitive 
cell contains two metal atoms which corresponds to 0.5 
ML coverage. The results for Al/Si(001) are elaborated 
here and the results for Ga and In are briefiy discussed. 
Two metal atoms in the cell may be placed in various 
ways on the surface. We consider several configurations 
to identify energetically the most favorable configuration. 
The configurations considered here are shown in the Fig. 
1. 

The metallic dimers configurations on reconstructed 
Si(OOl) shown in Figs. 1(a), 1(b) and 1(c) were consid- 
ered by Batra^-'^ . He found that the configuration in Fig. 
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FIG. 1: Top layer of reconstructed Si(OOl) surface. The 2x2 
cells closed by thick borders are shown for different configura- 
tions of metal atoms at 0.5 ML coverage. The unfilled circles 
represent the top layer Si atoms and the filled circles represent 
the metal atoms. 



1(b) was the most stable one followed by the configura- 
tion in Fig. 1(c). The configurations in Figs. 1(c) and 
1(f) were considered by Northrup et al2S and they con- 
cluded that configuration in 1(f) was energetically the 
most favorable one. The arrangement in 1(c) is referred 
to as the 'Ortho' orientation of metal dimers and that 
in 1(f) is the 'Para' orientation of metal dimers as de- 
scribed earlier by the Northrup et al^^. In addition we 
consider two more configurations and they are shown in 
Figs. 1(d) and 1(e) respectively. The relative energies 
and the distance between the neighbor atoms for all the 
configurations are given in table I. All the energies are 
measured with respect to the reference energy which at 
0.5 ML coverage of metal atoms correspond to the config- 
uration shown in Fig. 2(a). For calculating the reference 
energy, all the Si atoms are held fixed and the metal 
atoms are allowed to relax along the [001] direction. The 
reference energies for the cases of Al, Ga, and In are - 
131.625 eV, -130.714 eV, and -130.255 eV. respectively. 
It is clear from table I that the configuration shown in 
1(f) i.e., Al dimer placed between the Si dimer rows with 
'Para' orientation has the lowest energy and hence the 
most favorable configuration. Our results support the 
calculations first made by Northrup et and later by 
Brocks et al2^. We also observe that the next most favor- 
able configuration is that shown in 1(e) which also has 
'Para' orientation of Al dimer but the dimer is placed 
within a Si dimer row. However, the energy difference 
between these two configurations with 'Para' orientation 
of metal dimers is ~ 1 eV. This is because the dimer 
within the Si row [1(e)] has to reside at much higher po- 
sition compared to that placed between the dimer rows 
[1(f)] and hence has a weaker binding. We also note that 
the energy difference between the 'Para' and 'Ortho' ori- 
entation of dimers is even larger, the 'Para' being more 
stable '-^1.3 eV. 

To calculate the potential energy variation of the sys- 
tem as a function of the Al dimer orientation, we place 
the Al dimer between the Si dimer rows at a fixed height 
and rotate about the [001] direction. The result is shown 
in Fig. 3. This shows two minima at 'Ortho' and 'Para' 
orientations respectively and there is a small barrier as 
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FIG. 2: The configurations (a) and (b) are considered to cal- 
culate the reference energies at 0.5 and 1 ML metal coverages 
respectively. The unfilled circles and the filled circles repre- 
sent the top layer Si atoms and the metal atoms respectively. 
Here all the Si atoms are held fixed and the metal atoms are 
allowed to relax along the (001) direction. 



one goes from 'Ortho' to 'Para' orientation. Therefore, 
the Al dimers take 'Para' orientation which is strongly 
more favorable compared to the 'Ortho' orientation. This 
result also agrees with the observation in a recent exper- 
iment by Sakama et We note that though the Si 
dimers are stretched by 0.26 A the reconstruction of the 
Si surface is not lifted at 0.5 ML of coverage. 

We have also done similar calculations for Ga and 
In adsorption on the reconstructed Si(OOl) surface and 
found that the 'Para' orientation of Ga and In dimers 
are also favored with relative energies -3.29 eV and -2.90 
eV respectively. The results for Ga and In also agree 
with recent experimental results obtained by Sakama et 
aliSiiS. The surface pattern at 0.5 ML coverage is 2x2. 
It is also clear that all the patterns observed at low tem- 
peratures and below 0.5 ML coverage of Al, Ga and In on 
the reconstructed Si(OOl) can be explained by the 'Para' 
dimers as proposed by NorthrupS^, but contrary to the 
earlier proposal by Batr 
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Ortho Para 

FIG. 3: The variation of energy of a 2x2 reconstructed cell 
as the Al dimer is rotated from 'Ortho' to 'Para' orientation. 
The Al dimer is placed between the Si dimer rows of the 
reconstructed Si(OOl) surface. 

One focus of Batra's work^-'^, however, was to find a sta- 
ble metallic configuration for Al. He reported a zigzag 
nano-structure shown in Fig. 1(b) which he labeled as 
Rx(Hl, H2), as being stable with a local minimum in en- 
ergy. Our present calculations support this finding ex- 
cept that we get Ax = ±0.6 A in contrast to Batra's 
value of ±0.8 A. This overlayer structure of Al atoms is 
reminiscent of a uniform bond (CH)x structure with the 
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important difference that there is no direct Al-Al bond. 
The nearest neighbor Al-Al distance (~ 4 A) in the over- 
layer is much longer than the bulk Al-Al distance (~ 2.8 
A) and they "see" each other only through the delocal- 
ized electrons in the px-orbitals. In (CH)x, C-C atoms 
are not only coupled through 7r-electrons but they also 
form strong uniform cr-bonds leading to a semi-metallic 
structure. This structure becomes semi-conducting upon 
undergoing a bond alternation Peierls like transition to 
Double (D), Single (S), DS, DS ... repeating pattern of 
bonds. It supports the solitonSS, formation by creating a 
defective structure of the type .. DS, DS, SD, SD .. by 
having two single bonds adjacent to each other. The Al 
overlayer structure^ has two electrons of each Al atom 
tied to a surface Si atom in an AIX2 like fashion while 
the third electron of Al is in a free electron like state 
along the nanowire (y-direction). Here X represents a 
surface Si atom which nominally has a single electron 
in a dangling p^-orbital available for bonding with Al. 
Interatomic distance between Al and surface Si atoms, 
d(Al-X) 2.4 A, implies a strong bond. The delocalized 
Pz electron (one per Al atom) creates a one dimensional 
semi-metallic system. Naturally, under such situation 
metallic wires may be realized and it will be conduct- 
ing under small bias conditions. But then the standard 
Peierls distortion may come into play depending upon 
the strength of the surface bonds. Within the error of 
our calculations (estimated at ~ 7 meV) wc did not find 
any lowering of the total energy upon dimerization, un- 
like (Ci/)x. In fact we noted that a dimerization of Al 
atoms in the zigzag structure by an amount Ay — ±0.15 
A raised the total energy by ^ 0.1 eV. Thus the sur- 
face bonds formed by supported nanowires can lock the 
metallic structure in place. There is no possibility of a 
soliton formation in the Al overlayer structure because of 
the large Al-Al distance. However, it will be interesting 
to look for soliton behavior in other overlayers. 



B. 1 ML of coverage 

We consider the adsorption of four metal atoms on the 
Si(OOl) surface of a 2x2 unit cell which corresponds to 
1 ML coverage. Results for Al, Ga, and In on Si(GOl) 
are presented here. Various possible configurations con- 
sidered are shown in Fig. 4. Here calculations are done 
for both the ideal and the reconstructed Si surfaces to 
find energetically the most favorable configuration. The 
2x2 unit cells are bordered by thick lines in Fig. 4. The 
configurations considereed include Para-Para and Ortho- 
Ortho metal dimers on reconstructed (Figs. 4(a) - 4(d)) 
and Ideal (Figs. 4(e) - 4(h)) on Si(OOl) surface. In addi- 
tion, we also consider mixed dimers (Ortho and Para) on 
reconstructed surfaces (Figs. 4(i) - 4(1)). Mixed dimers 
on ideal Si surface were not stable and readily caused the 
silicon surface to reconstruct. 

The relative energies for the relaxed structures of Al, 
Ga and In on the reconstructed surfaces are given in table 
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FIG. 4: Metal atoms on top layer of reconstructed and ideal 
Si(OOl) surface. The 2x2 cells closed by thick borders are 
shown for different configurations metal atoms at 1 ML cov- 
erage. The unfilled circles represent the top layer Si atoms 
and the filled atoms represent the metal atoms (Al or Ga or 
In). 



II for four possible configurations and the same on the 
ideal surfaces are shown in table III. The energies for the 
configurations with mixed dimers are given in table IV. 
The reference energies at 1 ML coverage of metal atoms 
are obtained for the configuration shown in Fig. 2(b) 
where all Si atoms are held fixed the metal atoms are 
allowed to relax along the (001) direction. The reference 
energies for the cases of Al, Ga and In are -137.963 eV, 
-136.485 eV and -136.044 eV respectively. 

A general conclusion is that with increasing coverage 
beyond 0.5 ML, 'Ortho' dimers become part of the low 
energy configurations leading towards a 'Para' to 'Ortho' 
transition at 1 ML coverage. For In, mixed staggered 
dimers ('Ortho' and 'Para') give the lowest energy config- 
uration. For Ga, mixed dimers arc non-staggered, while 
for Al 'Para' to 'Ortho' transition of dimers is complete. 
This is an important conclusion because at or below 0.5 
ML 'Ortho' dimers were not favorable. Also, the puzzle 
that 'Ortho' and 'Para' dimers are both seen is resolved 
due to the preponderance of 'Ortho' dimers at increasing 
coverages. 

Let us first consider the structure consisting of Al 
dimers with 'Para' orientation as it was the most favor- 
able structure at 0.5 ML coverage. Two configurations 
are possible with 'Para' orientation of Al dimers: (i) the 
configurations shown in Figs. 4(a) and 4(e) on the re- 
constructed and ideal surfaces respectively and (ii) the 
configurations shown in Figs. 4(b) and 4(f) on the recon- 
structed and ideal surfaces respectively. After complete 
relaxation, we note that the final low energy structure is 
independent of the initial surface pattern (reconstructed 
/ideal). When we start from an ideal surface, the Si 
atoms dimerize and the low energy structure is the one 
with reconstructed Si surface. Same final structure is 
achieved when we start form the reconstructed surface. 
Also note that the relaxed system with the staggered 



5 



dimers configuration in 4(b) consisting of parallel dimers 
on the reconstructed surface is less favorable compared 
to that on the same reconstructed surface with the con- 
figuration in 4(a). Comparing 4(b) and 4(f) we find that 
the final energies of the systems differ by a large amount: 
staggered dimers on the ideal surface are least favorable. 

Next we consider the configurations in 4(c) and 4(d). 
These structures are obtained by just rotating the Al 
dimers in 4(a) and 4(b) by 90 degree {i.e., dimers now 
take 'Ortho' orientation). From table II it is clear 
that these configurations are more favorable compared 
to those in 4(a) and 4(b). We find that after complete 
relaxation, the total energies and the structures of these 
two configurations become independent of the initial sur- 
face pattern. 

In addition, we consider the configurations consisting 
of mixed dimers ('Ortho' and 'Para') as shown in Figs. 
4(i) - 4(1). Similar configurations on ideal surface are un- 
stable, the Si surface reconstructs itself and eventually 
leads to the configurations 4(i), 4(j) and 4(k). Again for 
Al, the configurations 4(j) and 4(1) are unstable. From ta- 
ble IV we see that the configuration with mixed Al dimers 
shown in Fig. 4(i) has a deep local minumum. This con- 
figuration is more favorable compared to the configura- 
tions consisting of only 'Para' orientation of Al dimers 
(i.e., Figs. 4(a) and 4(f)). 

However, comparing the energies for configurations 
4(i) and 4(g) we find that the lowest energy structure 
is the one shown in Fig. 4(g) where the reconstruction 
is completely lifted. All the Si bonds are saturated here. 
Since the reconstruction is completely lifted at 1 ML cov- 
erage, we conclude that the breaking of Si dimer bonds is 
initiated just beyond 0.5 ML of coverage. Most interest- 
ingly we note that an orientational transition for the Al 
dimers takes place as one goes from 0.5 ML coverage to 
1 ML coverage: The Al dimers change their orientation 
from 'Para' to 'Ortho' after breaking the Si dimer bonds. 
This is also supported by the charge density plot for the 
lowest energy configurations at 0.5 and 1 ML coverage as 
shown in Figs. 5(a) and 5(b) respectively. 
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FIG. 5: The charge density plot for: (a) the lowest energy 
structure at 0.5 ML coverage, (b) the lowest energy structure 
at 1 ML coverage. The charge density plane coincides with 
the Al layer. 



Earlier, it was argued that at 0.5 ML coverage all Si 
bonds are saturated and therefore beyond 0.5 ML cover- 
age Al atoms may reside on top of each other and form 3D 
island. One experimental study observed the 3D island 
formation beyond 0.5 ML of Al coverage^^. Our calcula- 
tion, however, suggests that for Al an orientational tran- 
sition can take place at 1 ML coverage and the surface 
structure becomes regular having a 1 x 2 pattern where all 
bonds of Si and Al are saturated due to removal of the 
reconstruction of the Si surface. This leads us to spec- 
ulate that further experiments on Al may seek out 1x2 
pattern at 1 ML coverage. Furthermore, at coverages be- 
tween 0.5 ML and 1 ML both the 'Para' and 'Ortho' Al 
dimers may coexist. This also clarifies why experiments 
at times have reported both 'Para' & 'Ortho' dimers be- 
cause precise coverage is a difficult parameter to quantify. 



(a) 



(b) 



Si — , „ _ . s 


i s 

i s 


Si Ga c 


S 



it: 


A -- -- Si 

i 


4l F 


1 . — ^i^^^j 



FIG. 6: The charge density plot for: (a) the lowest energy 
structure at at 1 ML coverage of Ga on Si(OOl), (b) the lowest 
energy structure at at 1 ML coverage of In on Si(OOl) 

Now we consider the Ga adsorption at 1 ML coverage. 
If we start with the reconstructed surface we find that 
the configuration in Fig. 4(i) with mixed dimers leads to 
the lowest energy structure where the reconstruction of 
the Si surface is not lifted during the relaxation process. 
The next best configuration is shown in Fig. 4(g) and 
consists of only 'Ortho' dimers. The structure remains 
ideal even after complete ionic relaxation. The configu- 
rations in Figs. 4(j) and 4(1) are unstable. Comparing 
final energies for various configurations on ideal and re- 
constructed surfaces we observe that the configuration 
in 4(i) is energetically the most favorable structure. Ex- 
periments indicate observing the 1x2 structure shown in 
4(g)ii. Here the reconstruction of the Si surface is lifted, 
Ga dimers prefer 'Ortho' orientation leading to a 1x2 
patten. Our calculations suggest that mixed dimers (Fig. 
4(i)) and 'Ortho' dimers (Fig. 4(g)) as two contenders for 
the energy favored structures. 

For the case of In we find that energetically the most 
favorable configuration is given in Fig. 4(j). This con- 
figuration corresponds to a 2x2 surface pattern. The 
configurations in Figs. 4(i) and 4(1) are unstable. More 
significantly, we note that the reconstruction of the sur- 
face is not lifted. We also find that the configurations 
consisting of 'Para' dimers (i.e.. Figs. 4(a) and a(b)) are 
degenerate and energetically next favorable structures. 
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The configurations in Figs. 4(a) and 4(b) correspond to 
the 2x1 and 2x2 surface structures respectively. Exper- 
imental results indicate the existence of both the 2x2 
and 2x1 surface structures near 1 ML coverage of In. 
The most and the next favorable structures from our cal- 
culations are consistent with the experimentally observed 
2x2 and 2x1 surface structures respectively. 

To understand the bonding of Ga and In atoms on 
Si(OOl), the charge density distributions are plotted on 
the plane of metal atoms for the most favorable configu- 
rations at IML coverage of Ga and In and they are shown 
in Figs. 6(a) and 6(b) respectively. We observe that to 
a first approximation metal atoms display strong intra- 
plane bondings amongst themselves leaving the substrate 
in reconstructed state. The charge distribution in Fig. 
6(b) also helps us to understand why the configuration 
4(i) is unstable for In. In being larger in size, the In atoms 
in 'Ortho' orientation of configuration 4(i) are pushed 
away from each other to have a suitable coupling among 
In atoms and thus favoring the configuration 4(j) with 
(mixed) staggered dimers. 

It is worth mentioning that one usually ends up with 
semi-metallic or non-metallic nanowires when supported 
on substrates as we saw above for the group III metals 
on Si(OOl). Very recently, there have been several im- 
portant studiesSSiSSi^SiSii^^ of monovalent Au chains on 
various high indexed Si surfaces dealing with mctallicity, 
Peierls distortion, charge density instability, fractional 
band fillings, and possibility of a spin charge separation 
due to the Luttinger-liquid phase. In fact, Himpsel and 
his group showed that Au nanowire supported on Si(553) 
is stable and it shows metallic behavioriSSiS. 



IV. SUMMARY 



We have done extensive total energy calculations for 
Al, Ga and In adsorption on the Si(OOl) surface at 0.5 and 
1 ML coverages. The 'Para' orientation of dimers is the 
most stable orientation at 0.5 ML coverage for Al, Ga and 
In, which supports the results of Northrup et al and also 
agrees with the recent experimental observations. All the 
patterns observed in various experiments at low temper- 
ature and below 0.5 ML coverages can be explained by 
the 'Para' orientation of metal dimers. An interesting 
zigzag semi-metallic structure similar to (CH)x is formed 
at 0.5 ML but is not the lowest energy structure. Soliton 
formation in overlayers is a fundamental topic for further 
study. 

At 1 ML coverage a tendency towards the 'Ortho' ori- 
entation of dimers is noted. Thus we predict an orienta- 
tional transition from 'Para' to 'Ortho' as we go from 0.5 
ML to 1 ML coverage for Al and possibly Ga on the Si 
surface. In addition, the reconstruction of the Si surface 
is completely lifted at 1 ML coverage of Al. These results 
lead us to conclude that the surface patterns at 0.5 ML 
coverage can be explained in terms of 'Para' orientation 
of metal dimers and those at intermediate coverages be- 
tween 0.5 and 1 ML may be explained by the co-existence 
of 'Para' and 'Ortho' orientation of dimers. For In mixed 
'Para' and 'Ortho' dimers give the lowest energy configu- 
ration that exhibit 2x2 surface pattern and the transition 
is not complete. 
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TABLE I: Relative energies of the relcixed 2x2 cell at 0.5 ML coverage of Al on the reconstructed Si(OOl) surface for various 
configurations. Various bond lengths are also given. The reference energy for this case is -131.625 eV corresponding to the 
configuration in Fig. 2(a). 



Configuration in figure 1 


Energy (eV) 


d(Si-Si) 


d(Si-Al) 


d(Al-Al) 


height 


1(a) 


-1.892 


2.52 


2.67 


3.84 


1.36 


1(b) 


-1.951 


2.46 


2.50 


4.01 


1.40 


1(c) 


-2.216 


2.78 


2.62 


2.64 


0.73 


1(d) 


-2.153 


2.59 


2.39 


2.60 


1.92 


1(e) 


-2.499 


2.45 


2.47 


2.62 


1.55 


1(f) 


-3.482 


2.44 


2.48 


2.67 


1.10 



TABLE II: Relative energies of a relaxed 2x2 cell at IML coverage of Al, Ga and In on a reconstructed Si(OOl) surface 
for four different configurations with either 'Ortho' or 'Para' dimers. "R" and "I" implies that the Si(OOl) surface becomes 

"Reconstructed" and "Ideal" respectively after complete ionic relaxation. The reference energies for the case of Al, Ga and In 
arc -137.963 cV, -136.485 cV and -136.044 eV respectively corresponding to the configuration in Fig. 2(b). 



Configuration in figure 4 


Al: Energy (eV) 


Ga: Energy (eV) 


In: Energy (eV) 


4(a) 


-2.680 (R) 


-2.528 (R) 


-2.077 (R) 


4(b) 


-2.371 (R) 


-2.448 (R) 


-2.077 (R) 


4(c) 


-3.092 (I) 


-2.148 (R) 


-1.441 (R) 


4(d) 


-2.868 (I) 


-2.187 (R) 


-1.441 (R) 



TABLE III: Relative energies of a relaxed 2x2 cell at IML coverage of Al, Ga and In on an ideal Si(OOl) surface for four different 
configurations with cither 'Ortho' or 'Para' dimers. "R" and "I" implies that the Si(OOl) surface becomes "Reconstructed" 
and "Ideal" respectively after complete ionic relaxation. The reference energies for the case of Al, Ga and In are -137.963 eV, 
-136.485 eV and -136.044 cV respectively corresponding to the configuration in Fig. 2(b). 



Configuration in figure 4 


Ah Energy (eV) 


Ga: Energy (eV) 


In: Energy (eV) 


4(e) 


-2.680 (R) 


-2.528 (R) 


-1.840 (R) 


4(f) 


-0.339 (I) 


-0.427 (I) 


0.099 (I) 


4(g) 


-3.092 (I) 


-2.739 (I) 


-0.389 (I) 


4(h) 


-2.868 (I) 


-2.631 (I) 


-0.364 (I) 



TABLE IV: Relative energies of a relaxed 2x2 cell at IML coverage of Al, Ga and In on a reconstructed Si(OOl) surface 
for four different configurations consisting of mixed ('Ortho' and 'Para') dimers. "R" implies that the Si(OOl) surface becomes 
"Reconstructed" after complete ionic relaxation. The configurations marked 'Unstable' transforms to one of other stable 
structure upon ionic relaxation. The reference energies for the case of Al, Ga and In are -137.963 eV, -136.485 eV and -136.044 

eV respectively corresponding to the configuration in Fig. 2(b). 



Configuration in figure 4 


Ah Energy (eV) 


Ga: Energy (eV) 


In: Energy (eV) 


4(i) 


-2.937 (R) 


-2.995 (R) 


Unstable 


4(j) 


Unstable 


Unstable 


-2.496 (R) 


4(k) 


-2.667 (R) 


-2.195 (R) 


-1.886 (R) 


4(1) 


Unstable 


Unstable 


Unstable 



